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1.  GOAL  OF  THE  PROJECT. 

Most  materials  are  either  crystalline  or  amorphous. 
Crystalline  materials  have  a  periodic  atomic  structure  while 
amorphous  materials  are  aperiodic.  Quasicrystals,  on  the  other 
hand,  have  an  atomic  structure  characterized  by  two  or  more 
periodicities  which  are  incommensurate  to  each  other.  Thus  the 
structure  has  long  range  orientational  and  positional  order,  but  no 
translational  symmetry.  This  departure  from  the  requirement  of 
translational  symmetry  allows  the  solids  to  have  symmetries  not 
possible  for  crystals,  such  as  icosahedral  symmetry  and  decagonal 
symmetry . 

New  symmetries  not  allowed  for  crystals  should  result  in  new 
symmetry  related  properties  so  far  unobserved  for  crystals.  The 
property  we  focussed  is  magnetism,  since  the  dependence  of  magnetic 
properties  on  the  crystalline  symmetry  is  well  known.  For 
instance,  we  have  shown  that  the  magnetic  anisotropy  energy  of  an 
icosahedral  solid  takes  the  form,  to  the  lowest  order, 


E(M)  -K6  [  Y%-  ( rpr  > 1/2  <  Yt~  *"<T5)  ] 


- 1  [231cos60-315cos40  +  lO5cos20-5+42cos0sin50cos5<J>] 
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where  are  the  spherical  harmonics  of  the  magnetization  direction 
( 6 ,  <p) .  Note  that  the  lowest  order  is  the  6-th  order,  with  4-th 
and  2nd  order  missing.  Since  the  magnitude  of  the  anisotropy 
energy  decreases  very  quickly  with  the  increasing  order,  the 
magnetic  anisotropy  energy  of  a  ferromagnetic  icosahedral 
quasicrystal  is  much  smaller  than  the  cubic  crystals  which  are 
dominated  by  the  4-th  order  terms,  and  hexagonal,  tetragonal  or 
orthorhombic  crystals  which  have  large  two-fold  terms.  Thus  the 
icosahedral  quasicrystalline  ferromagnets  should  exhibit  very  small 
magnetic  anisotropy  and  consequently  high  magnetic  permeability, 
and  may  provide  some  use  as  a  soft  magnetic  material.  The  goal  of 
the  present  project  has  been  to  develop  such  ferromagnetic 
quasicrystals  and  study  their  structure  and  properties. 


2.  ACHIEVEMENTS. 


At  the  time  we  submitted  the  research  proposal  of  this  project 
no  ferromagnetic  quasicrystal  has  been  known.  Together  with  Dr. 
O'Handley  of  MIT  and  his  colleagues  we  found  an  indication  that  a 
melt-quenched  (rapidly  cooled)  Co-Er  alloy  may  be  a  ferromagnetic 
quasicrystal.  Further  research,  howevei ,  has  shown  that  the  alloy 
was  merely  crystalline,  and  not  quasicrystalline.  We  made  further 
considerable  effort  by  sputter  depositing  the  Co-Er  and  Co-Er-Si 
alloy  system,  but  the  results  have  been  negative. 
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In  the  meantime  two  overseas  reports  claimed  an  observation  of 
ferromagnetism  in  quasicrystalline  solids.  The  first  was  by 
Chinese  researchers  on  Al-Fe-Ce  [1],  and  the  second  was  from  Japan 
on  Al-Gc-Mn  and  Al-Ge-Mn-Cu  [2].  The  magnetic  moment  reported  in 
Al-Fe-Ce  was  very  small,  suggesting  that  it  may  be  due  to  impurity 
phases.  At  the  same  time  amorphous  alloys  with  very  high  specific 
tensile  strength  were  produced  out  of  the  same  alloy  system  [3,4]. 
In  particular,  this  system  has  very  high  concentrations  of 
aluminum,  up  to  90  atomic  %.  Usually  amorphous  alloys  with  such  a 
high  concentration  of  the  primary  element  are  unstable  and  cannot 
be  produced  by  rapid  quenching,  since  they  often  violate  the  size 
rule  established  by  the  principal  investigator  [5].  It  was 
suspected  therefore  that  this  unusual  glass  forming  tendency  and 
the  formation  of  a  quasicrystal  may  have  something  in  common,  an  * 
we  decided  to  carry  out  the  neutron  and  x-ray  study  of  the  atomic 
structure  of  amorphous  Al-Fe-Ce.  As  expected  we  found  that  Fe 
strongly  modifies  the  local  structure  of  A1 ,  and  results  in  strong 
short  range  ordering,  as  discussed  below. 

The  other  alloy  system,  Al-Ge-Mn  and  Al-Ge-Mn-Cu,  shows  a 
larger  magnetic  moment,  and  apparently  is  more  stable.  However, 
its  magnetic  moment,  a,  is  quite  small  while  its  Curie  temperature, 
Tc,  is  relatively  high.  For  instance  for  Al40Ge25Mn25Cu10,  a  -  1.5 
emu/g  and  Tc  =  467  K.  Therefore  we  conjectured  that  while  it  was 
possible  that  this  magnetism  was  merely  due  to  a  small  amount  of 
the  second  phase,  possibly  crystalline,  such  as  the  grain  boundary 
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phase,  representing  uninteresting  possibility,  it  was  equally 
possible  that  the  system  was  actually  f errimagnetic  with  a  sizable 
local  moment,  and  the  smallness  of  the  total  moment  might  be  the 
consequence  of  partial  cancellation.  To  resolve  this  question  we 
studied  this  alloy  by  neutron  scattering.  The  result  shows  that 
the  second,  more  interesting  case  is  true,  and  we  indeed  made  the 
first  observation  of  the  magnetic  scattering  from  the  quasicrvstal . 
The  result  shows  that: 

1.  Strong  magnetic  order  can  exist  in  quasicrystals. 

2.  In  this  particular  case  the  magnetic  symmetry  is  lower  than 
the  lattice  (f.c.i.)  symmetry,  suggesting  that  it  is  either 
f errimagnetic  or  antiferromagnetic. 

3.  By  changing  the  antiferromagnetic  interaction  to  a 
ferromagnetic  one,  it  should  be  possible  to  produce  strongly 
magnetic  quasicrystals. 


2-1 .  Atomic  and  Magnetic  Structure  of  Ouasicrystalline 
Al40Ge25Mn2SCulo^ 


The  atomic  and  magnetic  structure  of  quasicrystalline 
Al40Ge„Mn25Cu10  was  studied  through  x-ray  and  neutron  scattering 
measurements.  The  melt-spun  flake  samples  were  provided  by  Prof. 
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Masumoto's  group  in  the  Institute  for  Materials  Research,  Tohoku 
University,  Sendai,  Japan.  The  x-ray  scattering  measurements  were 
made  using  synchrotron  radiation  at  the  X-7A  beamline  of  the 
National  Synchrotron  Light  Source  (NSLS)  of  Brookhaven  National 
Laboratory.  The  synchrotron  radiation  was  monochromated  by  a 
channel-cut  Si  crystal  at  25  keV,  and  the  scattered  x-rays  were 
detected  by  a  Ge  energy  sensitive  detector.  The  pulsed  neutron 
scattering  measurements  were  done  at  the  Special  Environment  Powder 
Diffractometer  of  the  Intense  Pulsed  Neutron  Source  (IPNS)  of 
Argonne  National  Laboratory,  and  magnetic  neutron  scattering 
measurements  at  the  reactor  of  the  National  Institute  for  Standard 
and  Technology  (NIST). 

The  indexing  of  the  diffraction  peaks  shows  that  there  exists 
an  extinction  rule,  that  the  six  indices  are  either  all  even  or  all 
odd,  similarly  to  the  extinction  rule  for  the  face  centered  cubic 
(f.c.c.)  structure.  Therefore,  this  solid  has  the  face  centered 
icosahedral  (f.c.i.)  structure,  just  as  quasicrystalline  Al6SCu2SFelt) 
[6].  We  have  shown  recently  that  the  structure  of  Al6SCu2SFe10  is 
given  by  the  projection  of  the  primitive  six-dimensional  hypercubic 
lattice  with  the  NaCl  type  atomic  decoration  using  the  projection 
window  inflated  by  r  compared  to  the  standard  window  to  produce  the 
3-d  Penrose  tiling  [7].  This  conclusion  was  obtained  by  comparing 
the  atomic  pair  distribution  function  (PDF)  determined  by  the 
synchrotron  x-ray  scattering  measurement  with  that  calculated  for 
the  model,  as  shown  in  Fig.  1.  This  is  the  first  example  of  the 
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primitive  lattice  quasicrystal,  since  we  also  found  that  other 
quasicrystals,  such  as  Al-Li-Cu  or  Al-Mn-Si,  have  a  structure 
projected  from  the  non-primitive,  decorated  6-d  lattice,  with  a 
high  density  of  local  phason  defects,  or  the  lattice  sites 
unoccupied  by  atoms.  We  suggest  that  the  primitive  nature  of  the 
6-d  lattice  is  the  origin  of  the  absence  of  phason  defects  in 
Al65Cu25Fe10 .  The  phason  defects  are  the  local  deviations  from  the 
perfect  quasicrystalline  structure,  and  linear  phason  defects  which 
broaden  the  diffraction  peaks  have  been  observed  in  all  the 

i 

icosahedral  solids  except  for  this  solid.  The  PDF  determined  for 
Al-Ge-Mn-Cu  indicates  that  the  same  structure  describes  this  solid 
as  well,  if  A1  is  replaced  by  random  mixture  of  A1  and  Ge,  Cu 
replaced  by  Mn,  and  Fe  by  Cu. 

i 

Non-spin-polarized  neutron  scattering  experiments  were 
performed  at  the  NIST  in  Gaithersburg,  MD  using  neutrons  with  2.6 
A  in  wavelength.  The  powder  sample  weighing  about  50  g  was 
contained  in  an  A1  tube,  and  placed  in  an  A1  low  temperature 

i 

i 

furnace.  We  found  that  the  height  of  the  diffraction  peaks  change 
ver>  significantly  by  annealing.  In  the  same  temperature  range 
Tsai  et  al.  [2]  found  a  large  exo-thermic  peak  in  the  differential 
scanning  calorimetry,  and  concluded,  without  other  evidences,  that 
it  was  due  to  crystallization.  However,  since  the  peak  positions 
do  not  change,  we  do  not  think  that  crystallization  is  taking 
place.  Instead,  what  is  happening  seems  to  be  the  increase  in  the 
6-d  NaCl  type  atomic  chemical  ordering,  and  the  modelling  certainly 
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confirms  this. 


Most  of  the  measurements  were  carried  out  on  the  sample 
annealed  at  275  C  for  30  min.  to  maximize  the  ordering.  The 
results  showed  that  the  many  diffraction  peaks  depend  upon 
temperature  much  more  than  expected  from  the  Debye-Waller  factor  as 
shown  in  Fig.  2.  In  particular  one  peak  at  Q  =  2.53  A"1  disappears 
above  270  K,  at  the  Curie  temperature  as  shown  in  Fig.  3. 
Therefore  this  peak  must  be  a  purely  magnetic  peak.  Thus  we  have 
for  the  first  time  observed  the  magnetic  neutron  scattering  from  a 
quasicrystal . 

The  extra  peak  at  2.53  A"1  can  be  indexed  as  (110000),  with  the 
quasicrystalline  lattice  constant  of  8.90  A,  thus  violating  the 
f.c.i.  extinction  rule.  Therefore  the  magnetic  sublattice  must  be 
breaking  the  f.c.i.  symmetry.  This  suggests  that  the  magnetic 
order  is  either  antiferromagnetic  or  f errimagnetic .  While 
ferromagnetic  sublattice  is  commensurate  with  the  lattice  and  does 
not  break  the  lattice  symmetry,  antiferromagnetic  or  f errimagnetic 
order  tends  to  double  the  lattice  periodicity.  We  are  planning  to 
".arry  out  spin-polarized  neutron  scattering  experiments  at  trie 
High-Flux  Reactor  of  Brookhaven  National  Laboratory  (under  the  NSF 
support)  when  the  reactor  resumes  operation  most  likely  early  next 
year,  to  separate  the  magnetic  peaks  more  clearly  from  the  nuclear 
peaks  and  to  determine  the  magnetic  structure.  It  is  not  clear  at 
this  moment  whether  the  system  is  truly  a  ferrimagnet,  or  a 
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partially  frustrated  antiferromagnet .  In  the  latter  case  a  perfect 
quasicrystal  would  have  no  magnetic  moment,  and  the  moment  is  due 
only  to  the  defects.  This  observation  suggests  that  if  by  chemical 
modification  the  anti ferromagnetic  exchange  interaction  could  be 
changed  to  ferromagnetic  interaction  there  is  a  good  possibility 
for  a  strong  ferromagnetic  quasicrystal. 

2 . 2  Atomic  Structure  of  Amorphous  Al-Fe-Ce. 

The  atomic  structure  of  melt-spun  ribbons  of  Al90FexCe10_x  was 
studied  by  pulsed  neutron  and  x-ray  scattering  measurements.  The 
samples  were  provided  by  Prof.  Poon's  group  at  the  University  of 
Virginia.  As  detailed  in  two  publications  on  this  subject  attached 
as  Appendices,  the  neutron  and  x-ray  scattering  measurements 
clearly  showed  that  the  Al-Fe  distance  in  this  alloy  system  is 
anomalously  short  compared  to  the  sum  of  the  atomic  radii,  about 
2.5  A  vs.  2.69  A.  This  is  most  likely  the  consequence  of  the 
strong  s-d  hybridization  between  A1  and  Fe.  This  strong  Fe-Al 
interaction  results  in  strong  compositional  as  well  as  topological 
short  range  ordering  around  Fe,  and  Fe  has  only  six  close  A1 
neighbors  at  2.5  A  and  4  more  slightly  distant  A1  neighbors  at 
about  2.8  A. 

Also  the  ordering  among  Fe  atoms  is  very  far  from  random.  The 
structure  factor  of  this  alloy  system  has  a  strong  pre-peak  at 
around  1.4  -  1.5  A"1,  which  extrapolates  to  about  1.6  A"1  for 
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Al<>0Fe10,  while  for  a  random  distribution  of  Fe  atoms  with  a  short 
range  ordering  due  to  the  repulsion  among  Fe  we  expect  to  have  a 
peak  at  1.3  A"1.  This  peak  position  is  close  to  that  of  a 
quasicrystalline  Al-Fe  [8],  suggesting  that  the  Fe  ordering  may  be 
close  to  that  in  quasicrystals.  Conversely  this  result  suggests 
that  the  strong  interaction  between  Fe  and  A1  (in  general 
transition  metals  and  A1 )  may  be  the  driving  force  behind  the 
formation  of  quasicrystals  of  this  class. 

On  the  other  hand  the  interaction  between  A1  and  Ce  appears  to 
be  quite  normal.  The  Al-Ce  distance  is  close  to  the  sum  of  the 
atomic  radii,  and  the  short  range  order  is  basically  that  of  the 
random  close  packing  ( DRP ) ,  with  the  repulsion  between  Ce  atoms. 
Thus  the  short  range  order  among  Ce  atoms  can  be  described  by  the 
dilute  DRP  model  of  Ce  atoms  in  the  amorphous  A1  matrix.  The  role 
of  Ce  in  this  glass  therefore  is  to  stabilize  the  glass  phase  by 
randomness,  following  the  mechanism  suggested  by  the  P.  I.  [5], 
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Figure  Captions: 

Fig.  1.  Atomic  pair  distribution  function  (PDF)  of  icosahedral 
Al65Cu25Fe10  determined  by  x-ray  scattering  (dashed  line) 
and  calculated  from  the  quasicrystalline  model  (solid 
line) . 

Fig.  2.  Neutron  diffraction  intensity  of  icosahedral 
Al40Ge25Mn25Cu10  at  10  K  (solid  line)  and  at  503  K  (above 
the  Curie  temperature,  dashed  line). 

Fig.  3.  Neutron  scattering  intensity  of  the  magnetic  peak  at  2.53 
A"1  at  various  temperatures.  Tc  is  at  194  C. 
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The  atomic  structure  of  liquid-quenched  amorphous  AlsoFeXeto-x  (x  =  5,7)  was  studied 
by  pulsed  neutron  and  x-ray  scattering.  The  atomic  pair-density  function  determined  by 
pulsed  neutron  diffraction  indicates  that  a  significant  portion  of  Al-Fe  distances  is 
anomalously  short,  while  some  part  of  the  Al-Al  distances  is  anomalously  long.  Both 
neutron  and  x-ray  scattering  showed  the  presence  of  a  prepeak  in  the  structure  factor. 
These  results  suggest  that  a  strong  interaction  between  A1  and  Fe  modifies  the  structure 
of  this  glass,  leading  to  chemical  and  topological  short-range  ordering. 


I  INTRODUCTION 

A  group  of  ductile  metallic  glasses  with  aluminum  up 
to  90  at.  %  were  recently  discovered  and  were  found  to 
have  excellent  strength  per  weight,  or  specific  strength.1,2 
Examples  of  these  glasses  are  Al-TM-RE  (TM  =  Fe, 
Co,  Ni,  etc.;  RE  =  Ce,  Y,  Gd,  etc.),  such  as  A^NisYg, 
Al9oFe5Ce5,  and  A^Fe^Gdo.1"3  The  tensile  strength 
of  AIwFe5Ce5  reaches  as  high  as  940  MPa,2,3  while  its 
physical  density  is  only  3.21  g/cm3.  Usually  metallic 
glasses  with  such  a  high  concentration  of  the  primary 
element,  in  this  case  Al,  cannot  be  formed  by  rapid 
quenching  as  they  most  often  fail  to  satisfy  the  atomic 
size  criterion  for  glass  formability,4  as  we  discuss  later. 
The  purpose  of  this  work  is  to  determine  the  origin  of 
the  glass  formability  through  the  study  of  their  atomic 
structure  using  pulsed  neutrons  and  synchrotron  x-rays. 
We  show  that  a  strong  interaction  between  Fe  and  Al 
leads  to  chemical  and  topological  short-range  ordering 
around  Fe,  and  this  contributes  to  the  ease  of  forma¬ 
tion  of  this  glassy  alloy. 

II.  EXPERIMENTAL  RESULTS 

Amorphous  ribbons  of  Al9oFe,Ceio-x(*  =  5,7) 
were  prepared  by  melt-spinning  in  a  partial  helium  at¬ 
mosphere  using  a  copper  wheel  of  20  cm  in  diameter, 
with  a  typical  circumferential  velocity  of  40  m/s.  In  gen¬ 
eral,  the  dimensions  were  15  /im  thick,  1  to  2  mm  wide, 
and  up  to  several  meters  long.  The  ribbons  were  very 
flexible  and  could  easily  be  bent  in  half  without  fractur¬ 
ing.2  The  pulsed  neutron  scattering  experiments  were 
carried  out  at  the  Special  Environment  Powder  Diffrac¬ 
tion  (SEPD)  spectrometer  of  the  Intense  Pulsed  Neutron 
Source  (IPNS)  of  the  Argonne  National  Laboratory. 
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The  x-ray  scattering  measurements  were  performed  at 
the  X-7A  beam-line  of  the  National  Synchrotron  Light 
Source  (NSLS)  of  the  Brookhaven  National  Laboratory. 

The  chemical  and  geometrical  atomic  short-range 
order  in  amorphous  materials  can  be  determined  by  the 
neutron  and  x-ray  scattering  experiments.5  For  multi- 
component  amorphous  systems,  three  formalisms  (Faber 
and  Ziman,6  Ashcroft  and  Langreth,7  and  Bhatia  and 
Thomton8)have  been  proposed  to  discuss  the  composi- 
tionally  resolved  structure  factors  to  describe  the  chemi¬ 
cal  and  geometrical  order.  In  the  following,  we  use  the 
Faber-Ziman  formalism.  The  Faber-Ziman  total  struc¬ 
ture  factor  in  the  case  of  neutron  scattering  is  defined  as: 

K*(Q)  «  (b)2S™(Q)  +  (b2)  -  (b)2  (1) 

where  /?*  is  the  coherent  scattering  intensity  per  atom, 
(b)  is  the  coherent  scattering  length  averaged  over  con¬ 
stituent  atoms,  (b2)  is  the  square  of  the  coherent  scatter¬ 
ing  length  averaged  over  constituent  atoms,  and  Q  = 
47 r  sin  0/A,  where  0  is  the  scattering  angle  and  A  is  the 
wavelength  of  the  incident  neutrons.  /«**  is  determined 
from  the  neutron  scattering  data  corrected  for  absorp¬ 
tion,  multiple  scattering,  and  Placzek  effect.910 

The  Faber-Ziman  total  structure  factor  for  an 
n-component  amorphous  system  is  expressed  as  the 
weighted  sum  of  the  partial  structure  factor  S0(Q)  asso¬ 
ciated  with  the  i-j  pair  correlation,  where  i(j)  is  the  in¬ 
dex  of  z-th(/-th)  components,  — -- 

sFZ(0  =  2  2.Cc;^s,(Q)  (2) 

i-l  ;-l  SP) 

where  c,  and  b,  are  concentration  fraction  and  coherent 
neutron  scattering  length  of  the  i-th  element,  and  n  is  the 
number  of  the  elements.  Note  that  St,  -  S,,.  The  partial 
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pair  density  function  representing  the  i-j  pair  correla¬ 
tions  derived  by  the  Fourier  transform  of  Sij{Q)  is: 

P„(r)  =  p o  +  ~  J  Q[S„(Q)  -  1]  sin  Qr  dQ 

(3) 

The  total  pair  density  function  p(r),  which  is  the 
Fourier  transform  of  SF2(£?)>  is  the  weighted  sum  of  the 
partial  pair  density  functions. 

p(r)  =  P0  +  ~  Jo’  Q[S?Z(Q)  -  1]  sin  QRdQ 
“iici cMpt(r)  (4) 

i-i  (o) 

The  same  procedure  applies  to  x-ray  data,  but  x-ray 
data  have  to  be  corrected  in  addition  for  polarization 
and  Compton  scattering,  and  normalized  to  obtain 
/a°\n  Also,  the  coherent  neutron  scattering  length  b 
should  be  replaced  by  x-ray  atomic  form  factor  f(Q )  in 
all  equations. 

The  neutron  and  x-ray  scattering  measurements 
were  performed  on  Al<wFeICeto-x  samples  with  x  -  5, 7. 
Figure  1  shows  the  total  structure  factor  of  Al9oFe7Ce3 
and  Al9oFe5Ces  from  the  neutron  scattering  experi¬ 
ments.  The  total  structure  factor  shows  a  double- hump 
main  peak  at  about  3  A-1  and  a  small  prepeak  at  about 
1.5  A"1.  The  total  structure  factor  from  the  x-ray  scat¬ 
tering  experiment  shows  similar  features,  as  shown  in 
Fig.  2.  A  prepeak  was  observed  also  for  amorphous 
Al9oY10  and  AIwYgNis.12  Table  I  lists  the  position  and 
area  of  the  prepeak  of  the  AlwFesCes  and  Al9oFe7Ce3 
samples  from  neutron  and  x-ray  experiments.  The  posi¬ 
tion  of  the  prepeak  shifts  to  higher  Q  as  the  concentra¬ 
tion  of  Fe  increases.  We  may  extrapolate  the  prepeak 


FIG.  1.  Total  structure  factor,  S{Q),  of  amorphous  Al^FesCej 
(solid  line)  and  A)«oFe7Cc)  (dashed  line)  determined  by  the  pulsed 
neutron  scattering  experiment. 


L  Q  .Ci/A)  ^ 

FIG.  2.  Total. Iow-to-medium  Q  portion  of  the.^tructure  factor, 
$(2),  of  AlwFesCes  (solid  line)  and  Al«Fe7Cej  (dashed  line)  deter¬ 
mined  by  the  x-ray  scattering  experiment. 


TABLE  1.  First  prepeak  of  total  structure  factor. 


Composition 

Type  of  experiment 

Peak-position 

Area 

AlqoFesCcj 

Neutron 

1.47 

0.099 

Al9oFc7Ccj 

Neutron 

1.53 

0.074 

AlwFcsCcj 

X-ray 

1.34 

0.170 

AlgoFciCej 

X-ray 

1.42 

0.132 

position  to  that  expected  for  the  compositions  Al^Feio 
and  Al9oCeI0  (i.e.,  AI^FezCejo-x  with  x  -  10,0).  The 
extrapolated  prepeak  positions  for  x  =  10  and  0  are 
1.58  ±  0.04  A"1  and  1.23  ±  0.09  A"1,  respectively.  The 
latter  agrees  with  the  prepeak  position  of  amorphous 
Al„Y,o,  1.3  A'1.12 

The  total  pair  density  function  (PDF)  of  AlgoFejCes 
and  Al9oFe7Ce3  from  neutron  experiments  is  shown  in 
Fig.  3.  The  first  peak  of  the  total  PDF  apparently  is  not 
a  single  peak,  but  consists  of  at  least  three  subpeaks. 
The  first  peak  of  the  radial  distribution  function  (RDF), 
47rr2p(r),  of  Al9oFe7Ce3  was  fitted  by  three  Gaussian 
functions,  as  shown  in  Fig.  4.  The  positions  and  areas 
for  the  three  subpeaks  are  listed  in  Table  II.  Also 
shown  in  Table  II  are  the  similar  data  for  the  first  three 
subpeaks  of  the  RDF  of  Al^FesCes  and  the  difference 
between  these  two  RDFs.  The  data  for  Al9oFe5Ce3  are 
derived  by  fitting  the  first  three  subpeaks  of  the  differ¬ 
ence  between  the  RDFs  of  Al9oFe7Ce3  and  Al^FesCes, 
as  shown  in  Fig.  5,  and  subtracting  them  from  the  data 
of  the  Al9oFe7Ce3  to  prevent  the  ambiguity  in  deter¬ 
mining  the  subpeaks.  Since  the  area  of  the  first  sub¬ 
peak  increases  with  the  Fe  concentration  and  Fe  is  the 
smallest  atom  in  this  alloy  system,  it  is  reasonable  to 
assign  the  first  subpeak  to  the  Al-Fe  correlation.  The 
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pair  density  function  representing  the  i-j  pair  correla¬ 
tions  derived  by  the  Fourier  transform  of  Sy(Q)  is: 

Pv(r)  =  po  +  rV  [  Q[Sti{Q)  -  1]  sin  Qr  dQ 

2 7TT  Jo 

(3) 

The  total  pair  density  function  p(r),  which  is  the 
Fourier  transform  of  SFZ(Q ),  is  the  weighted  sum  of  the 
partial  pair  density  functions. 

p(r)  =  Po  +  j  2[SFZ((2)  -  1]  sin  QRdQ 

=  2  20Cl~ip,l(r)  (4) 

The  same  procedure  applies  to  x-ray  data,  but  x-ray 
data  have  to  be  corrected  in  addition  for  polarization 
and  Compton  scattering,  and  normalized  to  obtain 
/*°Vl  Also,  the  coherent  neutron  scattering  length  b 
should  be  replaced  by  x-ray  atomic  form  factor  f{Q)  in 
all  equations. 

The  neutron  and  x-ray  scattering  measurements 
were  performed  on  AlooFe,Ce  io- r  samples  with  j:  =  5,7. 
Figure  1  shows  the  total  structure  factor  of  Al9oFe7Ce3 
and  AlgoFejCef  from  the  neutron  scattering  experi¬ 
ments.  The  total  structure  factor  shows  a  double- hump 
main  peak  at  about  3  A'1  and  a  small  prepeak  at  about 
1.5  A~l.  The  total  structure  factor  from  the  x-ray  scat¬ 
tering  experiment  shows  similar  features,  as  shown  in 
Fig.  2.  A  prepeak  was  observed  also  for  amorphous 
Al9oYlo  and  Al87Y8Ni5.12  Table  I  lists  the  position  and 
area  of  the  prepeak  of  the  Al9oFe5Ce5  and  Al9oFe7Ce3 
samples  from  neutron  and  x-ray  experiments.  The  posi¬ 
tion  of  the  prepeak  shifts  to  higher  Q  as  the  concentra¬ 
tion  of  Fe  increases.  We  may  extrapolate  the  prepeak 
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FIG.  1.  Total  structure  factor,  S[Q),  of  amorphous  AlwFesCej 
(solid  line)  and  AlooFeiCej  (dashed  line)  determined  by  the  pulsed 
neutron  scattering  experiment. 
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FIG.  2.  'Total. low-to-medium  Q  portion  of  the.^tructurc  factor, 
5(G),  of  AlwFcsCes  (solid  line)  and  Al9oFe7Cej  (dashed  line)  deter¬ 
mined  by  the  x-ray  scattering  experiment. 


TABLE  I.  First  prepeak  of  total  structure  factor. 


Composition 

Type  of  experiment 

Peak-position 

Area 

AlwFesCcs 

Neutron 

1.47 

0.099 

AlwFe7Cej 

Neutron 

!  c3 

C.074 

AlooFejCe5 

X-ray 

1.34 

0.170 

AlwFejCej 

X-ray 

1.42 

0.132 

position  to  that  expected  for  the  compositions  Al^oFeio 
and  AlwCeto  (i.e.,  AJ^FexCew-x  with  x  =  10,0).  The 
extrapolated  prepeak  positions  for  x  =  10  and  0  are 
1.58  ±  0.04  A'1  and  1.23  ±  0.09  A”1,  respectively.  The 
latter  agrees  with  the  prepeak  position  of  amorphous 
AIwYto,  1.3  A".11 

The  total  pair  density  function  (PDF)  of  AJ9oFe5Ce5 
and  Alc>oFe7Ce3  from  neutron  experiments  is  shown  in 
Fig.  3.  The  first  peak  of  the  total  PDF  apparently  is  not 
a  single  peak,  but  consists  of  at  least  three  subpeaks. 
The  first  peak  of  the  radial  distribution  function  (RDF), 
47 rr2p(r),  of  AJ9oFe7Ce3  was  Fitted  by  three  Gaussian 
functions,  as  shown  in  Fig.  4.  The  positions  and  areas 
for  the  three  subpeaks  are  listed  in  Table  II.  Also 
shown  in  Table  II  are  the  similar  data  for  the  first  three 
subpeaks  of  the  RDF  of  Al9oFe5Ce5  and  the  difference 
between  these  two  RDFs.  The  data  for  Al9oFe5Ce5  are 
derived  by  fitting  the  first  three  subpeaks  of  the  differ¬ 
ence  between  the  RDFs  of  Al9oFe7Ce3  and  Al9oFe5CeSl 
as  shown  in  Fig.  5,  and  subtracting  them  from  the  data 
of  the  Al9oFe7Ce3  to  prevent  the  ambiguity  in  deter¬ 
mining  the  subpeaks.  Since  the  area  of  the  first  sub¬ 
peak  increases  with  the  Fe  concentration  and  Fe  is  the 
smallest  atom  in  this  alloy  system,  it  is  reasonable  to 
assign  the  first  subpeak  to  the  AI-Fe  corre  lukiV/al.  The 
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FIG.  5.  Total  pair  oensity  function  of  AlwFcjCcj  (solid  line)  and 
AI<»Fe7Cej  (dashed  line)  derived  from  the  neutron  scattering 
experiment. 


Al  coordination  around  Fe  calculated  from  the  area  of 
the  first  subpeak  is  6.2  and  6.3  for  Al^FesCes  and 
Al9oFe7Ce3,  respectively,  which  are  identical  within  the 


r  .(A) 

FIG.  4.  The  first  peak  of  the  radial  distribution  function  of 
AI«Fe7Cej  fitted  by  three  Gaussian  functions  (dashed  lines). 
Areas  and  positions  of  the  three  Gaussian  functions  are  listed  in 
Table  II. 


TABLE  II.  Positions  and  areas  of  first  three  subpeaks  of  the  RDF 
of  AlooFe7Cej  and  of  the  difference  between  the  RDFs  of 
AlooFe7Ce}  and  AlooFejCe*. 


Al<*)Fe7Cej  Al«>oFc7Cej— AlgoFejCcs 

AlwFesCe: 

Area* 

Position 

Area 

Position 

Area 

1st  subpeak 

2.49 

1.88 

2.51 

0.49 

1J9 

2nd  subpeak 

2.81 

6.21 

2.72 

0.29 

5.92 

3rd  subpeak 

3.25 

2.93 

3.30 

-0.63 

3.56 

r*  _  j. 

CL 

'  r- 

•  «  W  1 

*  ■  *  "V 

FIG.  5.  Difference  between  the  radial  distribution  function  of 
Al9oFe7Cej  and  the  radial  distribution  function  of  AlsoFcsCej,  fit¬ 
ted  by  three  Gaussian  functions  (dashed  lines).  Areas  ana  posi¬ 
tions  of  the  three  Gaussian  functions  are  listed  in  Table  II. 

margin  of  error,  justifying  the  assignment  of  this  sub¬ 
peak  to  the  Fe-Al  correlation.  The  position  of  this  sub¬ 
peak  (2.49  A),  however,  is  less  than  the  sum  of  the 
atomic  radii  of  AI  (1.43  A)  and  Fe  (1.26  A)  by  about 
0.2  A.  The  peak  positions  of  the  second  subpeak 
(2.81  A)  and  the  third  subpeak  (3.25  A)  are  close  to  the 
expected  Al-Al  distance  (2.86  A)  and  Al-Ce  distance 
(3.25  A),  respectively. 

The  average  coordination  number  of  an  impurity 
A  atom  embedded  in  the  glass  of  B  element  can  be  es¬ 
timated  accurately,  within  the  dense  random  packing 
(DRP)  model,  using  the  principle  of  two-dimensional 
packing  of  the  atomic  surface.13  The  average  coordi¬ 
nation  number  is  a  function  of  the  size  ratio  of  these 
two  atoms: 

Mab  =  V2d(4(Wvab) 

=  4 77(1  -  V5/2)/{l  -  [rA(rA  +  2 r0)]1/2/ 

(rA  +  rB)}  (5) 

where  rj2 d  is  the  surface  packing  ratio,  rA  and  rB  are  the 
atomic  radii  of  the  elements  A  and  B,  respectively,  and 
c oab  is  the  solid  angle  to  view  B  atom  from  the  central 
A  atom.  The  computer  simulation  showed  that  this  pre¬ 
diction  is  quite  accurate  over  a  wide  range  of  size  ratios.14 
We  will  refer  the  coordination  number  NAB  calculated 
from  Hq.  (5)  as  the  value  predicted  by  the  DRP  model. 
For  the  Al-Ce  correlation  such  a  calculation  predicts  a 
contribution  of  1.08  to  the  area  of  the  third  peak  for 
Al9oFe7Ce3,  which  is  much  less  compared  to  the  experi¬ 
mental  value  of  2.93.  Therefore,  there  must  be  contribu¬ 
tions  from  other  correlations.  The  most  probable  source 
is  the  Al  atoms  in  the  vicinity  of  Fe  having  long  Al-Al 
distances.  As  we  discussed  earlier,  the  Fe  atom  has 
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about  6  close  AI  atom  neighbors.  If  these  AI  atoms  are 
evenly  packed  around  Fe,  the  distance  between  these 
AI  atoms  should  be  much  larger  than  twice  the  radius 
of  the  AI  atom,  in  the  neighborhood  of  3.3  A,  thus  con¬ 
tributing  to  the  third  subpeak.  If  we  assume  the  re¬ 
maining  contribution  of  the  area  of  the  third  subpeak 
ctimes  from  the  elongated  Al-Al  distances,  we  can  cal¬ 
culate  the  AI  coordination  around  AI  with  long  (3.3  A) 
distances  to  be  2.66,  which  is  quite  reasonable  for 
6  contacting  AI  atoms.  A  similar  estimate  made  for  the 
third  peak  of  the  Al9oFe5Ce5  indicates  the  AI  coordina¬ 
tion  around  AI  of  elongated  distance  is  2.27,  which  is 
slightly  smaller  than  that  of  Al#)Fe7Ce3,  but  this  may 
be  due  to  the  underestimation  of  the  area  of  the  third 
subpeak  which  is  overlapping  with  the  next  peak. 

Possible  origins  of  the  second  subpeak  of  RDF  are 
the  Al-Al  correlation  and  the  second  nearest  AI  neigh¬ 
bor  around  Fe.  The  partial  pair  density  function  pai,aiW 
within  the  1st  nearest  neighbor  subshell  will  decrease 
when  introducing  Fe  and  Ce.  Fe  and  Ce  atoms  occupy 
part  of  the  nearest  neighbor  subshell  of  AI  atoms  and 
decrease  the  number  of  AI  atoms  contacting  this  AI 
atom.  Fe  also  elongates  the  distance  of  AI  atoms  con¬ 
tacting  it  and  moves  some  of  the  nearest  neighbor  cor¬ 
relation  away  from  the  first  subshell  at  2.8  A.  We  can 
quantitatively  evaluate  the  contributions  by  summing 
all  the  effects.  For  simplicity,  we  examined  the  dif¬ 
ference  of  the  second  subpeak  of  the  RDF  of  Al^FesCes 
and  Al9oFe7Ce3.  There  is  2%  increase  in  the  Fe  num¬ 
ber  concentration  and  2%  decrease  in  the  Ce  number 
concentration  when  changing  from  AUoFesCes  to 
Al9oFe7Ce3.  The  calculation  showed  that  the  effects  of 
introducing  Fe  and  Ce  are  almost  equal  and  cancel 
each  other.  Therefore  the  increase  in  the  second  sub¬ 
peak  when  changing  from  Al9oFe5Ce5  to  Al9oFe7Ce3  is 
most  probably  due  to  the  second  nearest  neighbor  of  AI 
around  Fe.  With  this  assumption  the  calculated  second 
nearest  neighbor  coordination  of  AI  around  Fe  at  a  dis¬ 
tance  of  approximately  2.8  A  is  about  3.8.  Thus  there 
are  roughly  10  AI  atoms  surrounding  Fe  within  3  A, 
and  about  6  of  them  are  contacting  Fe,  just  as  in  the 
case  of  crystalline  compound  AI3Fe.15 

III.  DISCUSSION 

A  large  number  of  alloy  systems  were  found  to  form 
metallic  glasses  when  they  were  rapidly  quenched  from 
the  liquid  state.16  The  question  of  glass  formability— in 
other  words,  which  factors  would  influence  the  ease  of 
glass  formation— has  been  discussed  extensively,. and  a 
number  of  factors  have  been  proposed  to  be  the  glass 
formability  factors.4,17,18  They  include  the  atomic  size 
ratio,4,17  heat  of  alloy  formation,17  and  valence  electron 
concentration.18  Among  them,  the  atomic  size  ratio4  ap¬ 
pears  to  be  the  one  with  a  least  number  of  exceptions.19 
In  that  sense  it  could  be  identified  as  the  most  impor¬ 


tant  glass  formability  factor.  The  differential  atomic 
volume  ratio  was  found  to  be  inversely  proportional  to 
the  minimum  solute  concentration  necessary  to  form  a 
stable  glass  phase  by  rapid  quenching.4  The  approxi¬ 
mate  relationship  is  given  by 


Crain; 
B 


0.1 


(6) 


where  C£in  is  the  minimum  atomic  concentration,  of 
solute  atoms  B  in  the  matrix  of  atoms  A  needed  to  form 
a  stable  glass,  and  VA  and  VB  are  the  atomic  volume  of 
matrix  A  and  solute  B.  For  A  =  AI  and  B  =  Ce,  CB  in 
is  0.1;  while  for  A  =  AI  and  B  =  Fe,  CBm  is  0.35.  The 
size  difference  creates  atomic  level  stresses  which  cause 
the  topological  instability  of  a  crystalline  solid  solution 
and  facilitates  glass  formation.4  This  equation  could  be 
generalized  to  ternary  systems  by  summing  the  prod¬ 
ucts  of  the  minimum  atomic  concentration  and  differ¬ 
ential  atomic  volume  ratio  of  each  solute  element: 


—  0.1 


(7) 


where  B  and  C  are  two  solute  elements  and  A  is  the 
matrix  atoms.  This  equation  was  first  used  successfully 
to  describe  the  glass-forming  range  in  the  Co-TM-B 
system  (TM  =  Ti,Zr,Hf,Cr,Mo,W,V,Nb,Ta,Mn).20 
The  calculated  values  of  the  sum  of  the  products  of 
minimum  atomic  concentration  and  differential  volume 
ratio  for  Al9oFe5Ce5  and  Al9oFe7Ce3  are  0.067  and 
0.052,  respectively.  These  values  are  much  less  than  the 
expected  value,  0.1,  indicating  that  the  concentrations 
of  Fe  and  Ce  are  much  less  than  required  by  the  equa¬ 
tion.  This  is  mainly  due  to  the  small  atomic  size  differ¬ 
ence  between  AI  and  Fe. 

The  shortening  of  the  Al-Fe  distance  indicates  a 
strong  interaction  between  AI  and  Fe  which  might  be 
the  cause  of  high  glass  formability.  The  small  number 
of  AI  coordination  around  Fe  and  the  anomalously 
short  Fe-Al  distance  suggests  that  the  effective  Fe  ra¬ 
dius  in  AI  may  be  much  smaller  than  the  regular  metal¬ 
lic  radius,  due  to  tnis  interaction.  We  may  then  use 
Eq.  (5)  to  define  an  effective  size  ratio  according  to  the 
coordination  number  of  AI  around  Fe  (6.3).  From  the 
effective  size  ratio  and  the  experimentally  measured 
Al-Fe  nearest  neighbor  distance,  we  can  calculate  the 
effective  radius  of  Fe  and  AI  and  estimate  the  Al-Al 
nearest  neighbor  distance.  The  effective  rA/rB  calcu¬ 
lated  by  Eq.  (5)  is  0.47,  which  implies  that  the  effective 
radii  of  Fe  and  AI  are  0.80  A  and  1.69  A,  respectively. 
Thus  the  distance  between  the  adjacent  AI  atoms 
around  the  Fe  is  3.38  A,  which  is  quite  close  to  the  third 
subpeak  position  of  3.25  A.  The  effective  size  ratio 
also  predicts  the  values  of  the  sum  of  the  products  of 
minimum  atomic  concentration  and  differential  volume 
ratio  of  the  Al9oFe5Ce5  and  Al9oFe7Ce3  systems  to  be 
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0.097  and  0.094,  which  is  again  quite  close  to  0.1.  Thus 
this  assumption  of  the  reduced  effective  size  of  Fe  is 
internally  consistent,  and  explains  very  well  why  the 
unusually  small  concentration  of  the  solute  (Fe)  atoms 
can  stabilize  the  glass.  The  strong  Al-Fe  interaction 
and  the  short  Al-Fe  distance  suggested  here  are  consis¬ 
tent  with  the  large  negative  deviation  from  volume  ad¬ 
ditivity  of  the  Al-Fe  alloy  in  the  liquid  or  solid  state.21 
It  is  of  interest  to  explore  in  the  future  whether  the 
deviation  from  volume  additivity  is  generally  accompa¬ 
nied  by  unusual  interatomic  distances  or  not. 

This  rather  drastically  reduced  effective  size  of  Fe 
and  the  increased  size  of  Al  may  be  explained  from  the 
electronic  point  of  view.  The  d-orbitals  of  Fe  and  the 
s-p  orbitals  of  AJ  should  hybridize  strongly  to  form  a  de¬ 
localized  d-band,  which  will  be  filled  by  the  transfer  of 
electrons  from  the  s-p  states  as  a  result  of  the  s-p  band 
being  pushed  up  in  energy  by  the  d-band,  just  as  in 
the  case  of  transition  metal-metalloid  alloys.22  Conse¬ 
quently,  the  Fe-Ai  distance  is  reduced  due  to  the  in¬ 
crease::  Fe-Al  bonding  strength,  while  the  delocalized 
d-electrons  increase  the  Al-Al  distance  around  Fe. 
Perhaps  reflecting  the  symmetry  of  the  d-orbitals,  the 
s-p  states  of  only  six  AI  atoms  hybridize  with  the  Fe 
d-states  and  shorten  the  Fe-Al  distance,  while  the  re¬ 
maining  four  Al  atoms  are  only  weakly  in  contact  with 
Fe  at  the  distance  of  2.8  A.  A  related  but  slightly  differ¬ 
ent  factor  also  contributes  to  this  effect.  Al  has  three 
valence  electrons  per  atom  and  has  a  high  electron  den¬ 
sity  with  a  high  Fermi  level.  Fe  metal  has  7.4  d-orbital 
electrons  whose  Fermi  level  is  lower  than  that  of  AI. 
The  difference  in  Fermi  energy  might  cause  a  large 
amount  of  transfer  of  electrons  from  AI  to  Fe,  but  in 
metallic  systems  such  a  large  amount  of  internal  polar¬ 
ization  is  unstable.  Instead,  the  effective  size  of  Fe  will 
be  reduced  to  increase  the  Fermi  level,  while  that  of  Al 
will  be  increased  to  balance  the  flow  of  charges. 

Another  point  which  merits  further  discussion  is 
the  extrapolated  position  of  the  prepeak  of  the  total 
structure  factor  of  AlooFexCeio-xC*  =  0,10).  In  a  binary 
alloy  system  in  which  the  minority  atoms  are  repelling 
each  other  but  more  or  less  uniformly  distributed,  there 
will  be  a  subpeak  in  the  pm{r)  at  a  position  2r0  where  B 
represents  a  minority  element  and  r0  is  the  average 
nearest  neighbor  distance  between  the  minority  atoms. 
The  value  of  r0  can  be  calculated  by  assuming  a  single¬ 
component  dense  random  packed  system  with  an  atomic 
radius  r0  whose  number  density  is  equal  to  the  number 
density  of  the  minority  atoms  in  the  binary  alloy  sys¬ 
tem.  This  subpeak  will  be  transformed  to  a  prepeak  in 
the  partial  structure  factor  5Bb(C?)  at  2.57 r/r0  due  to  the 
law  of  Fourier-sine  transform,  the  inverse  relation  of 
Eq.  (3).  The  subpeak  in  the  pBB(^")  will  emerge  as  a 
small  prepeak  in  the  total  structure  factor  S(Q )  because 
of  the  small  concentration  of  the  minority  atoms  and 


thus  the  long  distance  between  them.  For  Al^Feio  and 
Al90Cei0,  the  prepeak  position  is  calculated  to  be 
1.28  A*1  and  1.22  A-1,  respectively.  The  extrapolated 
position  of  the  prepeak  of  AlwCeio  is  1.23  ±  0.09  A-1, 
while  the  prepeak  position  in  AI90Y10  is  1.3  A'1;  thus 
both  of  these  fit  the  value  predicted  by  the  dense  ran¬ 
dom  packing  model  quite  nicely.  This  suggests  that  the 
local  structure  of  amorphous  Al-Ce(Y)  might  be 
described  well  by  the  dense  random  packing  model, 
provided  that  the  Ce-Ce  potential  is  repulsive.  For 
Al9oFe!0,  however,  the  extrapolated  position  of  the  pre¬ 
peak  (1.58  ±  0.04  A-1)  of  the  structure  factor  is  quite 
different  from  the  value  predicted  by  the  DRP  model 
(1.28  A-1).  On  the  other  hand,  the  position  of  the  pre- 
peak  is  in  the  close  vicinity  of  one  of  the  prominent 
diffraction  peaks  of  quasicrystalline  Al-Fe  ([110001]  re¬ 
flection  at  about  1.6  A-1).23  This  suggests  that  the  local 
structure  of  AI  around  the  Fe  might  be  similar  to  the 
quasicrystalline  structure,  and  Fe-Al  clusters  may  be 
packed  with  some  icosahedral  short-range  order.  In  fact, 
a  quasicrystalline  icosahedral  phase  has  been  reported 
in  rapidly  cooled  Al-Fe-Ce.24  Conversely,  it  is  possible 
to  speculate  that  a  strong  interaction  such  as  the  one  ob¬ 
served  here  may  be  necessary  to  stabilize  the  quasicrys- 
talline  structure. 

IV.  CONCLUSION 

The  atomic  structure  of  Al^FeXeio-x  (x  =  5,7) 
was  studied  by  neutron  and  x-ray  scattering  measure¬ 
ments.  Analysis  of  the  radial  distribution  functions 
shows  that  the  Fe  atoms  are  surrounded  by  about  10  AI 
atoms  within  3  A,  with  6  of  them  being  in  close  contact 
with  Fe.  The  distance  between  the  Al-Fe  pair  in  close 
contact  (2.49  A)  is  less  than  the  sum  of  the  atomic 
radii,  suggesting  a  strong  interaction  between  Al  and 
Fe.  Also,  some  of  the  Al-Al  distance  is  elongated.  The 
glass  formability  with  such  a  high  concentration  of  Al 
and  the  elongated  Al-Al  distances  can  be  phenomeno¬ 
logically  explained  by  introducing  an  effective  atomic 
size  ratio  (0.47)  between  Fe  and  Al,  which  is  signifi¬ 
cantly  reduced  from  the  ratio  of  the  metallic  radii  of  each 
element.  The  drastic  departure  of  the  effective  atomic 
size  of  Fe  in  the  matrix  of  Al  from  the  regular  metallic 
radius  is  most  likely  due  to  the  consequence  of  the  strong 
s-d  hybridization.  The  analysis  of  the  prepeak  of  the 
total  structure  factor  indicates  that  Ce  atoms  are  basi¬ 
cally  randomly  distributed  in  the  Al  glass  except  that 
they  are  repelling  each  other,  while  a  strong  composi¬ 
tional  and  geometrical  order  exists  around  Fe  which 
may  be  similar  to  the  local  order  in  quasicrystals. 
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ANALYSIS  OF  THE  RESULTS.  The  total  pair  density  function  ^(r)  which  is  the  Fourier 
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